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Abstract The deflection of a ray of light passing close to a gravitational 
mass, is generally calculated from the null geodesic which the light ray ( 
photon) follows. However, there is an alternate approach, where the effect 
of gravitation on the ray of light is estimated by considering the ray to be 
passing through a material medium. Calculations have been done in this 
paper, following the later approach, to estimate the amount of deflection 
due to a static non-rotating mass. The refractive index of such a material 
medium, has been calculated in a more rigorous manner in the present work 
and the final expression for the amount of deflection calculated here is claimed 
to be more accurate than all other expressions derived so far using material 
medium approach. Based on this expression, the amount of deflection for a 
sun grazing ray has been also calculated. 

1 Introduction 

The gravitational deflection of light is one of the important predictions of the 
General Theory of Relativity (GTR) proposed by Einstein, which plays a key 
role in understanding problems related to Astronomy, Cosmology, Gravita- 
tional Physics and other related branches. 

Newtons theory of universal gravitation had already predicted that the 
path of any material particle moving at a finite speed is affected by the pull 
of gravity. By the late 18th century, it was possible to apply Newtons law 
to compute the deflection of light by gravity. Cavendish commented briefly 
on the gravitational deflection of light in the late 1700s and Soldner gave a 
detailed derivation in 1801. 
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The idea of bending of hght was revived by Einstein in 191 1 and the quan- 
titative prediction for the amount of deflection of light passing near a large 
mass (M) was identical to the old Newtonian prediction, d — 2GM/{c^rQ), 
where tq is the closet distance of approach and in this case approximately 
the solar radius. It wasn't until late in 1915, as Einstein completed the 
general theory, he realized his earlier prediction was incorrect and the angu- 
lar deflection should actually be twice the size he predicted in 1911. This 
was subsequently confirmed by Eddington in 1919 through an experiment 
performed during the solar eclipse. 

The exact amount of deflection for a ray of light passing close to a gravi- 
tational mass can be worked out from the null geodesic, which a ray of light 
follows [1,2,3]. 

The deflection of a light ray passing close to a gravitational mass can be 
alternately calculated by following an approach, where the effect of gravita- 
tion on the light ray is estimated by considering the light ray to be passing 
through a material medium with a value of refractive index decided by the 
value of gravitational fleld [4] . 

The concept of this equivalent material medium was discussed by Balazs 
[5] as early as in 1958, to calculate the effect of a rotating body, on the 
polarization of an electromagnetic wave passing close to it. Plebanski [6] 
had also utilized this concept in 1960, to study the scattering of a plane 
electromagnetic wave by gravitational fleld, where the author mentioned that 
this concept of equivalent material medium was first pointed out by Tamm 
[7] in 1924. A general procedure for utilizing this concept, for defiection 
calculation has been worked out by Felice [8]. Later this concept was also 
used by Mashoon [9,10], to calculate the deflection and polarization due to the 
Schwarzschild and Kerr black holes. Fischbach and Freeman [11], derived the 
effective refractive index of the material medium and calculated the second 
order contribution to the gravitational deflection. In a similar way Sereno 
[12] has used this idea, for gravitational Icnsing calculation by drawing the 
trajectory of the ray by Fermat's principle. More recently Ye and Lin [13], 
emphasized the simplicity of this approach and calculated the gravitational 
time delay and the effect of lensing. 

On the other hand, the calculation of higher order deflection terms, due to 
Schwarzschild Black hole, from the null geodesic, has been performed recently 
by Iyer and Fetters ( [14] and references their in). Using null geodesies, 
gravitational lensing calculations have been done by a number of authors in 
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past [15,16]. 

With the above background, in the present work, we follow the mate- 
rial medium approach, to calculate a more accurate value for the deflection 
term due to a non- rotating sphere ( Schwarzschild geometry ). It is claimed 
that the present calculated value will be more accurate than all other values 
calculated in past, using material medium concept. 

2 The effective refractive index and the tra- 
jectory of Hght ray 

As discussed earlier, the gravitational field infiuences the propagation of elec- 
tromagnetic radiation by imparting to the space an effective index of refrac- 
tion n(r)[4]. 

For a static and spherically symmetric gravitational field, the solution of 
Einstein's Field Equation was given by K. Schwarzschild in 1961, which is as 
follows [4]: 

ds^ = (1 - '^ydt'' - r\sin''edcj>^ + de'') - -^Y^T^ (1) 

where = called Schwarzschild Radius, which completely defines the 
gravitational field in vacuum produced by any centrally-symmetric distribu- 
tion of masses. The above equation can be expressed in an isotropic form by 
introducing a new radius co-ordinate ( p ) with the following transformation 
equation [4] 

p=l[(r-^)+rVV-r,)V^] (2) 
OR 

r = P(l + 'ff (3) 
4p 

The resulting isotropic form of Schwarzschild equation will be now: 

= - (1 + 'f)\dp^ + p\sinHdct>' + dO')) (4) 

l + rg/{Ap) 4p 
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Now in spherical co-ordinate system the quantity {dp^ + p'^{sin'^9d(tP' + 
dO"^)) has the dimension of square of infinitesimal length vector d~^ . 

By setting ds — 0, the velocity of hght can be identified from the ex- 
pression of the form ds^ — f{p)dt^ — d~^^, as v{p) — f{p)- Therefore the 
velocity of light in the present case ( characterized by Schwarzschild radius 
Tg ) can be expressed as : 

(1- r)c 

But this above expression of velocity of light is in the unit of length p per 
unit time. We therefore write 



v{r) = v{p) 



dp 

Substituting the value of p from Eqn (2) in Eqn.(6), we get: 



v{r) 



,Tj2-2p 
'rg/2 + 2p^ 

,r3/2-((r-^) + ri/2(r-r-^i/2 



9 J J^2^ 



r. 



,/2 + ((r-^)+rV2(^_^Ji/2) 



rg-r-r^'-[r-rgj ■ 
c(r - rg) 



(7) 



Therefore the refractive index n(r) at a point with spherical polar co- 
ordinate (r), can be expressed by the relation: 

n{r) = = ^— (8) 
v{r) r — rg 
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At this stage the entire problem, can become a problem of geometrical 
optics, where we have to find the trajectory of a light ray travelling in a 
medium, whose refractive index has spherical symmetry. The trajectory of 
the light ray and the center of mass ( source of gravitational potential) will 
together define a plane. The equation of such a ray in plane polar co-ordinate 
system {r,9) can be written as [17]: 



The trajectory is such that n{r).d always remains a constant, where d is 
the perpendicular distance between the trajectory of the light ray from the 
origin and the constant is taken here as A [17]. In our present problem the 
light is coming from infinity ( r = — oo) and it is approaching the gravitational 
mass, which is placed at the origin and characterized by Schwarzschild radius 
Tg. The closest distance of approach, for the approaching ray is h and the 
ray goes to r = oo, after undergoing certain amount of deflection (A0). 

Here, the parameter h can be replaced by solar radius tq. When the 
light ray passes through the closest distance of approach (ie r = 6 or ^q), 
the tangent to the trajectory becomes perpendicular to the vector 1^ ( which 
is "^q)- Therefore, we can write A — n{rQ)rQ. The trajectory of the hght 
ray had been already constructed before hke this, by Ye and Lin [13] and the 
value of deflection (A0), can be written as : 



However, Ye and Lin [13], had in our opinion used a value of refractive 
index n(r) which was approximated and somewhat ad hoc. Fischbach and 
Freeman [11] also in their attempt to calculate a more accurate value of 
deflection, considered terms only up to second order in the expression for 
refractive index. However, in our attempt to do so we shall avoid making 
any such approximation in the following. We denote the above integral in 
Eqn. (10) by / and write 






(10) 
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O ^\l^n{rQ).rQ 



dr 



n[rr 



= n{rQ)rQ ^ 



r^(n(r).r)2 - (n(rQ).rQ)^ 



o 



= n{rQ)rQ ^ 



dr 



O ^2 



1 Q 

(1-^)2 (1 ^) 



(11) 



Now we change the variable to x — ^ and introduce a quantity a — 
We also denote n{rQ) by hq. Accordingly we write: 



—X '^Vgdx 



a 



J a 



(a(l-a))2 

37 f 



xr^ 



riQTQ ra 

Jo 

f 

Jo 



{x{l-x)y (a(l-a))2 

dx 



'3 (a(l-a))2 



1 



1 



1 - 



(a(l-a))2 



(12) 



For our convenience we can denote the quantity l/(a(l — a)) by D. This 
also implies 



(13) 
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However, the Integral / can not be solved as a standard integral at this 
stage. We split the above Integral, as a sum of two Integrals and proceed as 
follows: 



O^Ou/"" {l-2x)dx ^ f'' xdx 
Jo Jl- D'^xHl-xY Jo .fl 



^1 - D^x^{l - x)2 Jo ^l-DH^{l-xf 



rg Jo Jl-D^x\l-xY ^9 Jo Ji-n^x^l-x 



\2 



= {^^)h + {^^)I2 (14) 

^9 ^9 

where Ii and I2 are used to denote the above two integrals. Now we can 
identify 

Tg 1 — a a a(l — a) 

Changing the variable from x to y = Dx{l — x), we can write D{1 — 2x)dx = 
dy. Accordingly the upper and lower limits x = and x = a change to y = 
and y = Da{l — = ^{i-a) '^^^ ~ ^ Therefore for the first part in 

Eqn (14) we can write : 



O^^O _ /■« D{l-2x)dx 



rg Jo ./l 



-D^x^{l-xf 
dy 



VI - y''^ 

= 7r/2 (15) 
Therefore, from Eqn (10), one may write the amount of defiection as: 



re 

A(/)= 2 / 

Jr,- 



„ ^ l( n[r).r p _ ^ 



— TT 
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^9 ^9 

= 7r + 2(-^=^)/2-7r 

^9 

^9 ^1 - L>2x2(l - xY 

Thus the gravitational bending for a ray of hght grazing the static grav- 
itational mass ( with Schwarzschild radius Vg ) with the closest distance of 
approach tq can be expressed as: 

A0^2L>/% ^"^^ (17) 
•^0 Jl- D^x^il-xY 



The above expression for gravitational deflection has been obtained from 
the Schwarzschild Equation ( Eqn(l)), with out applying any approximation 
at any stage. Owing to this, it is claimed that this expression of bending 
is more exact as compared to all other expressions derived till today, using 
equivalent material medium concept. However, the integration of the quantity 
in Eqn (17)), involves some complicated algebraic expressions containing 
Elliptical functions. Using Mathematica, we obtain the following expression 
after integration : 



/ 



xdx 



l-D^x'^il-x) 



[VD + y/D^)E - (2/D^)F 



(18) 



where E = E{p, q^) is the Elliptic Integral of first kind and F = F{—q,p, (f) 
is Incomplete Elliptic Integral of Third kind. The arguments p,?^,-q,p,Q'^ are 
expressed by the following mathematical relations: 



P 



arcsm 



{^fP^ - VdT4)(VD^+ {2x - l)^/D) 
\ {^/D^ + ^/D~+I)(^/D^ - (2x - 1)VD) 



(19) 
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(20) 

Finally we can write the expression for gravitational deflection (A0) of 
the hght ray, due to a static mass (r^) with the closest distance of approach 
tq as : 

^ ^ UVDW^E^^EMr (21) 

where the value of D is given by Eqn. (13) as D = , ^ — r and a — 

Tg/vQ. Eqn. (21) is a general expression for bending of light, where tq 
can be replaced by the closest distance of approach of the light ray. This 
mathematical expression for deflection, derived here is claimed to be more 
accurate than all other expressions derived so far using material medium 
approach and it is equally valid for strong field. For a Sun grazing ray, 
we can take the closest distance of approach as equal to solar radius which 
is tq = 695, 500 km and Schwarzschild radius corresponding to the mass 
of Sun as = 3 km. We, therefore, get a = {j-g/r-Q) = 1/231,833 and 
D = 231,834. Finally, we get a value of A0 = 8.62690^10"^ radians or 
1.77943 arc sec. This value of gravitational deflection suffered by a Sun 
grazing ray, is claimed to be more accurate than all other values obtained in 
past. 
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